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Based on a sample of rj mesons produced in the radiative decay J/if —► 7 r\ in 1.31 x 10 9 J/ip 
events collected with the BESIII detector, the decay rf —> u;e + e~ is observed for the first time, 
with a statistical significance of 8 <r. The branching fraction is measured to be Blrf —► Loe + e~) = 
(1.97 ± 0.34(stat) ± 0.17(syst)) x 10 —4 , which is in agreement with theoretical predictions. The 
branching fraction of rj —»• cwy is also measured to be (2.55±0.03(stat) ±0.16(syst)) x 10 2 , which is 
the most precise measurement to date, and the relative branching fraction ^ is determined 



to be (7.71 ± 1.34(stat) ± 0.54(syst)) x 10“ 3 . 
PACS numbers: 12.40.Vv, 14.40.Be, 13.20.Jf 
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I. INTRODUCTION 

The main decays of the v( meson [1] fall into two dis¬ 
tinct classes. The first class consists of hadronic decays 
into three pseudoscalar mesons, such as r / —> rjTrir, while 
the second class has radiative decays into vector particles 
with quantum number J PC = 1 , such as rj' —> 77,757, 

or wy. Model-dependent approaches for describing low 
energy mesonic interactions, such as vector meson dom¬ 
inance (VMD) [2] , and the applicability of chiral pertur¬ 
bation theory [2] can be tested in rf decays. 

It is of interest to study the decay 77' —> Ve + e~ (V 
represents vector meson) which proceeds via a two-body 
radiative decay into a vector meson and an off-shell pho¬ 
ton. The electron-positron invariant mass distribution 
provides information about the intrinsic structure of the 
r)’ meson and the momentum dependence of the tran¬ 
sition form factor. Recently, BESIII reported the mea¬ 
surement of 7 / -A TT + Tr~e + e~ [3], which is found to be 
dominated by 77' — >pe + e~, in agreement with theoretical 
predictions [2, 4]. 

Based on theoretical models [2, 5], the branching frac¬ 
tion of 77' —y ue + e~ is predicted to be around 2.0xl0 -4 , 
but until now there has been no measurement of this 
decay. A sample of 1.31 x 10 9 J/ij) events (2.25 x 10 8 
events [6] in 2009 and 1.09 x 10 9 [7] in 2012) has been 
collected with the BESIII detector and offers us a unique 
opportunity to investigate rj 1 decays via J/t/> —> 777'. In 
this paper, the observation of 77' -A ue + e“, the analysis 
of the decay 77' -A wy, and the ratio of their branching 
fractions are reported. 


II. DETECTOR AND MONTE CARLO 
SIMULATION 

The BESIII detector is a magnetic spectrome¬ 
ter located at the Beijing Electron Positron Collider 
(BEPCII, [8]), which is a double-ring e + e _ collider with 
a design peak luminosity of 10 33 cm -2 s _1 at a center-of- 
mass energy of 3.773 GeV. The cylindrical core of the 
BESIII detector consists of a helium-based multilayer 
drift chamber (MDC), a plastic scintillator time-of-flight 
system (TOF), and a Csl (Tl) electromagnetic calorime¬ 
ter (EMC), which are all enclosed in a superconducting 
solenoidal magnet providing a 1.0 T (0.9 T for the 2012 
run period) magnetic field. The solenoid is supported 
by an octagonal flux-return yoke with modules of resis¬ 
tive plate muon counters (MUC) interleaved with steel. 
The acceptance for charged particles and photons is 93% 
of the full 47 t solid angle. The momentum resolution 
for charged particles at 1 GeV/c is 0.5%, and the reso¬ 
lution of the ionization energy loss per unit path-length 


( dE/dx ) is 6%. The EMC measures photon energies with 
a resolution of 2.5% (5%) at 1 GeV in the barrel (end- 
caps). The time resolution for the TOF is 80 ps in the 
barrel and 110 ps in the end-caps. Information from the 
TOF and dE/dx is combined to perform particle identi¬ 
fication (PID). 

The estimation of backgrounds and the determinations 
of detection efficiencies are performed through Monte 
Carlo (MC) simulations. The BESIII detector is mod¬ 
eled with GEANT4 [9, 10]. The production of the J/if> 
resonance is implemented with the MC event generator 
KKMC [11, 12], while the decays are simulated with EVT- 
GEN [13]. Possible backgrounds are studied using a sam¬ 
ple of ‘inclusive’ J/tfj events of approximately the equiv¬ 
alent luminosity of data, in which the known decays of 
the J/ip are modeled with branching fractions being set 
to the world average values from the Particle Data Group 
(PDG) [1], while the remaining decays are generated with 
the LUNDCHARM model [14]. For this analysis, a sig¬ 
nal MC sample (6.0 x 10 5 events), based on the VMD 
model and chiral perturbation theory [2] for J/ip -A 777', 
77' -A coe + e~, cu — > 7r°7r + 7r _ , 7 r° -A 77 , is generated to op¬ 
timize the selection criteria and determine the detection 
efficiency. 


III. ANALYSIS OF J/ip -4 7rj' 

In this analysis, the r/ meson is produced in the radia¬ 
tive decay J/t/> —> 777'. The ui meson is observed in its 
dominant 7T + 7r _ 7r° decay mode, and the 7r° is detected 
in 7 r° -A 77 . Therefore, signal events are observed in 
the topology 77777 r + 7 r _ for the r/ -A wy mode, and 
7777r + 7r _ e + e _ for 77' —> we + eU We apply the follow¬ 
ing basic reconstruction and selection criteria to both 
channels: 

We select tracks in the MDC within the polar angle 
range | cos0| < 0.93 and require that the points of clos¬ 
est approach to the beam line be within ±20 cm of the 
interaction point in the beam direction and within 2 cm 
in the plane perpendicular to the beam. 

Photon candidates are reconstructed by clustering sig¬ 
nals in EMC crystals. At least four photon candidates 
are required, and the minimum energy of each must be 
at least 25 MeV for barrel showers (|cos0| < 0.80) and 
50 MeV for endcap showers (0.86 < | cos 6*| < 0.92). To 
exclude showers due to the bremsstrahlung of charged 
particles, the angle between the nearest charged track 
and the shower must be greater than 10°. To suppress 
electronics noise and energy deposits unrelated to the 
event, the EMC cluster time is restricted to be within a 
700 ns window near the event start time. 
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A. ri' -A w 7 

For the decay 77 ' -A ury, two particles with opposite 
charge are required. No particle identification (PID) is 
used, and the two tracks are taken to be positive and 
negative pions from the ui. 

A four-constraint (4C) kinematic fit imposing energy- 
momentum conservation is performed under the hypoth¬ 
esis of J/ip -A 77777 r + 7 r - . If there are more than four 
photons, the combination with the smallest xiy 77 
is retained. Events with y 2 , < 80 are retained 

/v 'y^'Y'77r ' 7r 

for further analysis. Since J/ip -A 777 ' is a two-body de¬ 
cay, the radiative photon carries a unique energy of 1.4 
GeV. Hence the photon with maximum energy is taken 
as the radiative photon, and its energy is required to be 
greater than 1.0 GeV. The photon pair combination with 
77 invariant mass closest to the 7r° mass is considered as 
the 7T° candidate in the final state, and its invariant mass 
must satisfy |M( 77 ) — M, r o| < 0.015 GeV/c 2 , where M n 0 
is the world average value of the 7 r° mass [1]. With these 
requirements, the decay 77 ' -a wy is observed in the dis¬ 
tribution of M(tt° 7t + 7t _ 7 ) versus M(7r°7r + 7r _ ), shown in 
Fig. 1. Besides the region of interest in Fig. 1 , there is 
a vertical band around the ui mass region, which comes 
from J/ip -a W 77 , urn 0 and ui7r°Tr° background, while a 
horizontal band also exists around the if mass region, 
which comes from J/ip -A 777 ', if -a r/n + Tr~ and 7 p°. 

To improve the mass resolution, as well as to better 
handle the background in the vertical band around the ui 
mass region and horizontal band around the if mass re¬ 
gion, we determine the signal yield from the distribution 
of the difference between M ( 7 r° 7 r + 7 r _ 7 ) and M (n°Tr + Tr~). 
The backgrounds in the vertical and horizontal bands do 
not peak in the signal region, which is demonstrated by 
the inclusive MG sample, as shown by the histogram in 
Fig. 2. 



M(7r°7r + ^') (GeV/c 2 ) 


FIG. 1. Distribution of M( 7 r° 7 r + 7 r 7 ) versus M(n 0 'K + 'K ) 
from data. 

To determine the signal yield, an unbinned maxi¬ 
mum likelihood fit to the mass difference M ( 7 r° 7 r + 7 r _ 7 ) — 


M( 7 r° 7 r + 7 r _ ) is performed, in which the signal shape is 
described by the MC shape convoluted with a Gaus¬ 
sian function to account for the difference in resolu¬ 
tion between data and MC simulation, and the back¬ 
ground is described by a 3rd-order Chebychev polyno¬ 
mial. 33187 ± 351 vf -a wy signal events are obtained 
from the fit, whose curve is shown in Fig. 2. With the 
detection efficiency, (21.87 ± 0.02)%, obtained from MC 
simulation, the branching fraction, (2.55 ± 0.03) x 10~~ 2 , 
listed in Table I, is determined. 



FIG. 2. Distribution of the mass difference M(n°n + n~'y) — 
M(jv°n + n~). The dots with error bars are data, the his¬ 
togram shows the MC simulation of inclusive J/ip decays. The 
solid curve represents the fit results, and the dashed curve is 
the background determined by the fit. 


B. rj' —> uje^e 

For rf -a we + e~ decay, candidate events with four 
well-reconstructed charged tracks and at least three pho¬ 
tons are selected. The charged track and good photon 
selections are exactly the same as described above. 

To select candidate events and select the best photon 
combination when additional photons are found in an 
event, the combination with the smallest xIc+pid * s re ~ 

tained. Here x\ c+vm = xlc + Ej=i XpidC?) is the sum 
of the chi-squares from the 4C kinematic fit and from 
PID, formed by combining TOF and dE/dx information 
of each charged track for each particle hypothesis (pion, 
electron, or muon). If the combination with the small¬ 
est xic+PiD corresponds to two oppositely charged pi¬ 
ons and an electron and positron, and has xlc < 80, the 
event is kept as a 7777 r + 7 r~e + e _ candidate. As in the 
analysis of rf —> wy, the selected photon with maximum 
energy is taken as the radiative photon, and its energy 
is required to be greater than 1.0 GeV. The other two 
photons are further required to be consistent with a 7r° 
candidate, |M(yy) — < 0.015 GeV/c 2 . 

With the above selection criteria, MC simulation shows 
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FIG. 3. (a) Distribution of the distance of the reconstructed e + e“ vertex from the 2 axis, R xy , where the dots with error bars 
are data, the solid histogram is signal MC simulation, and the dotted histogram is MC simulation of rf —>• 077 . (b) Distribution 
of M(e + e~) versus R xy , where the requirement of R xy < 2 cm is indicated as the vertical line, (c) Distribution of M(e + e~) 
with the requirement R xy < 2 cm, where the dots with error bars are data and the solid histogram is signal MC simulation. 


that background peaking under the signal comes from 
J/7/7 —> 777', 77' —>• W7, with the 7 from the 77' decay sub¬ 
sequently converting to an electron-positron pair. The 
distribution of the distance from the reconstructed ver¬ 
tex point of an electron-positron pair to the 2 axis, de¬ 
fined as R xy , is shown in Fig. 3 (a). As expected from 
MC simulation of J/ip —>- 777', 77' —>• 077, the peaks 
around R xy = 3 cm and R xy = 6 cm match the posi¬ 
tion of the beam pipe and the inner wall of the MDC, 
respectively, as shown in Fig. 3 (a). From the distri¬ 
bution of M(e + e~) versus R xy and the M(e + e~) pro¬ 
jections, shown in Figs. 3 (b) and (c), the requirement 
of R xy < 2 cm can cleanly discriminate signal from the 
background. The number of peaking background events 
from 77' — > 077 that still survive is estimated to be 2.6±0.3 
from MC simulation taking the branching fraction of 
J/7/7 — > 777', if — > 077 from this analysis, where the er¬ 
ror is statistical. This background will be subtracted in 
the calculation of the branching fraction of 77 ' —> uje + e~. 

With all the above selection criteria being applied, 
the scatter plot of M(7r°7r + 7r _ e + e _ ) versus 
is shown in Fig. 4 (a), where the cluster in the rj' 
and a; region corresponds to the decay 77 ' —»• 07e + e _ . 
The 77 ' and w peaks are clearly seen in the distribu¬ 
tions of M( 7 T° 7 r + 7 r _ e + e _ ) (Fig. 4 (b)) and M(n°TT + TT~) 
(Fig. 4 (c)), respectively. 

The same selection is applied to the inclusive MC sam¬ 
ple of 1.2 x 10 9 J/ 7/7 events to investigate possible back¬ 
ground channels. The corresponding normalized distri¬ 
butions of M(7r°7r + 7r _ e + e _ ) and M(n°TT + TT~) are shown 
as the histograms in Fig. 4 (b) and (c). One of the 
dominant backgrounds is from events with multiple 7r° 
in the final state with one 7 r° undergoing Dalitz decay to 
7e + e _ . Another important background, rj' —> n + n~r], 
77 — > TT°TT + n~ with the pion pair from the rj’ decay 
misidentified as an electron-positron pair, produces an 
accumulation at the low mass region in the distributions 
of M(7r°7r + 7r _ e + e _ ) and ), and at the high 

mass region in M(7r°7r + 7r _ e + e _ )-M(7r°7r + 7r~), which is 
shown as the shaded histograms in Fig. 4 (b), (c) and (d), 
normalized with the branching fraction from the PDG. 


The distribution of M( 7 r° 7 T + 7 r _ e + e _ ) - is 

shown in Fig. 4 (d). From this study of the inclusive MC 
sample, no peaking background events are expected. 

To determine the r( —> we + e“ yield, an unbinned max¬ 
imum likelihood fit to M(7r°7r + 7r _ e + e _ ) — 
shown in Fig. 5, is performed. The signal component is 
modeled by the MC simulated signal shape convoluted 
with a Gaussian function to account for the difference 
in the mass resolution between data and MC simula¬ 
tion. The shape of the dominant non-resonant back¬ 
ground 77' —>• is derived from the MC simulation, 

and its magnitude is fixed taking into account the de¬ 
cay branching fraction from the PDG [1]. The remaining 
background contributions are described with a 2 nd-order 
Chebychev polynomial. The fit shown in Fig. 5 yields 
66 ± 11 77' —> ue + e~ events with a statistical signifi¬ 
cance of 8 cr. The statistical significance is determined 
by the change of the log-likelihood value and of the num¬ 
ber of degrees of freedom in the fit with and without the 
77' —► U 7 e + e“ signal included. 

To determine the detection efficiency, we produce a sig¬ 
nal MC sample in which 77 ' —>• 07e + e“ is modeled as the 
decay amplitude in Ref. [2] based on the VMD model. Af¬ 
ter subtracting the peaking background events and tak¬ 
ing into account the detection efficiency of (5.45±0.03)%, 
the branching fraction of rj 1 —> uie + e~ is determined to 
be (1.97 ± 0.34) x 10 -4 . This is summarized in Table I. 


TABLE I. Signal yields, detection efficiencies and the branch¬ 
ing fractions of r\ —> 077 and 7/ — > uje + e~. The first errors 
are statistical, and the second are systematical. 


Decay mode 

Yield 

£ (%) 

Branching fraction 

77' —> 077 

77 ' —> u 7 e + e~ 

33187 ± 351 
66 ± 11 

21.87 

5.45 

(1.97 ±0.34 ±0.17) x 10 “ 4 
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FIG. 4. (a) Distribution of M(7r°7r + 7r _ e + e _ ) versus M(Tr°n + Tv~). (b) Invariant mass spectrum of 7r°7r + 7r _ e + e _ . (c) Invariant 

mass spectrum of 7r°7r + 7r _ . (d) Distribution of M(n°n + n~e + e~) — M(tt° n + n~). The solid histogram represents the remaining 
events from the inclusive MC sample, and the shaded histogram shows misidentified events from the background channel 
rf —¥ r/n + Tv~ normalized by using the branching fractions from the PDG [1]. 



FIG. 5. Distribution of M(7r°7r + 7 r _ e + e _ )— M(7r°7r + 7 r _ ) and 
the fit results. The crosses show the distribution of data. The 
dash-dotted line represents the rf —> Tv + n~g component, and 
the dotted curve shows the background except rf —> 7r + 7r ~ij. 


IV. SYSTEMATIC UNCERTAINTIES 

In this analysis, the systematic uncertainties on the 
branching fraction measurements mainly come from the 
following sources: 

a. MDC Tracking efficiency 

The tracking efficiencies of pions and electrons have 


been investigated using clean samples of J/if —> pn, 
if 1 -A Tr + n~J/ip 1 and J/if> —>• e + e~(jFSR)- Following the 
method described in Ref [15], we determine the difference 
in tracking efficiency between data and simulation as 1% 
for each charged pion and 1.2% for each electron. There¬ 
fore, 2% is taken as the systematic error of the tracking 
efficiency for rf —> wy with two charged tracks, and 4.4% 
for rf —>• we + e“ with four charged tracks. 

b. PID efficiency 

For rf —> u)e + e~, PID is used when we obtain xIc+p/d 
of every combination for each event. The decay J/ip —> 
7 r + 7r - 7T°, with n° —> ye+e - is used as a control sample to 
estimate the difference between data and MC with and 
without applying x^id t° identify the particle type. The 
difference, 3.8%, is taken as the systematic uncertainty 
from PID for the decay rf —>• we + eU 

c. Photon detection efficiency 

The photon detection efficiency has been studied in 
J/if -A prr decays in Ref. [15]. The difference between 
data and MC simulation is determined to be 1% per pho¬ 
ton. Therefore, 4% and 3% are taken as the systematic 
uncertainties, respectively, for the two analyzed rf de¬ 
cays. 

d. Kinematic fit 

The angular and momentum resolutions for charged 
tracks are significantly better in simulation than in data. 
This results in a narrower distribution in MC than 
in data and introduces a systematic bias in the efficiency 
estimation associated with the 4C kinematic fit. The 
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TABLE II. Summary of systematic uncertainties (in %) for 
the branching fraction measurements. 


Sources rf 

—»• oje + e 

77 ' — > uj’y 

B(w'— XjtJ'y) 

MDC tracking 

4.4 

2.0 

2.4 

Photon detection 

3.0 

4.0 

1.0 

PID 

3.8 

— 

3.8 

Kinematic fit 

1.8 

0.5 

1.9 

7 conversion subtraction 

1.0 

— 

1.0 

Background uncertainty 

3.7 

2.9 

4.7 

Form factor uncertainty 

1.3 

- 

1.3 

7T° mass window 

1.4 

1.4 

— 

J/f > total number 

0.8 

0.8 

— 

B{J/%p -*■ 777 ') 

3.1 

3.1 

- 

B(uj 7 r° 7 r + 7 r) 

0.8 

0.8 

- 

Total 

8.7 

6.4 

7.0 


difference can be reduced by correcting the track he¬ 
lix parameters of the simulated tracks, as described in 
detail in Ref. [16]. In this analysis, a clean sample of 
J/ip —> 7 T + 7 r _ 7 r°, 7 T° —> 7 e + e _ is selected to study the 
difference of the helix parameters of pions and electrons 
between data and MC simulation. The helix parame¬ 
ters of each charged track are corrected so that xlc fr° m 
MC simulation is in better agreement with that of data. 
With the same correction factors, the kinematic fit is per¬ 
formed for the signal MC events and the xlc is required 
to be less than 80. By comparing the numbers of selected 
signal events with and without the correction, we deter¬ 
mine the change in detection efficiencies to be 0.5% and 
1.8%. These are taken as the systematic uncertainties for 
rf —> uj'y and ?/ — > uie + e~ , respectively. 

e. 7 conversion event veto 

In the analysis of rj' —y uje + e~, the large contamina¬ 
tion of 7 conversion events from the decay r/' —y uj'y is 
effectively removed by the requirement of R xy < 2 cm. 
To estimate the uncertainty associated with this require¬ 
ment, we select a clean sample of J/ip —> 7 r + 7 r _ 7 r° with 
7 r° —> 7 e + e _ . The efficiency corrected signal yields with 
and without the R xy criterion differ by 1.0%, which is 
taken as the systematic uncertainty. 

/. Background 

The non-peaking background uncertainties in each 
channel are estimated by varying the fit range and chang¬ 
ing the background shape in the fit, and they are deter¬ 
mined to be 2.9% for 77 ' —> uj'y. To reduce the statisti¬ 
cal uncertainty for 77 ' — > ue + e“, we use the background 
shape from the inclusive MC sample, and the maximum 
change of the branching fraction, 3.6% is taken as the un¬ 
certainty from the non-peaking background. In order to 
evaluate the background uncertainty from rj' —> r]'ir + 'ir~ 
in the analysis of the r[ —y we + e _ decay, to, we perform 
an alternative fit by varying its contribution according to 
the uncertainty from branching fractions of J/ip —> 77 / 
and its cascade decays. We also vary the selection effi¬ 
ciency of this background channel as determined by the 
MC sample, and find that the total difference in the signal 
yield is about 0.3%, which can be ignored. In addition, 


the change in the number of peaking background events 
from 77 ' — y uj'y due to a difference of the 7 conversion 
ratio between MC and data leads to an uncertainty of 
1.0% on the branching fraction of 77 ' —> we + e~. The to¬ 
tal background uncertainties from these sources are listed 
in Table. II. 

g. Form factor 

The nominal signal MC model is based on the ampli¬ 
tude in Ref. [2] To evaluate the uncertainty due to the 
choice of the form factors in the determination of the 
detection efficiency, we also generate MC samples with 
other form factors in Ref. [2], e.g ., the monopole and 
dipole parameterizations. The maximum change of the 
detection efficiency, 1.3%, is regarded as the systematic 
uncertainty from this source. 

h. 7 r° mass window requirement 

The uncertainty from the 7 r° mass window requirement 
due to the difference in the mass resolution between data 
and simulation is estimated by comparing the difference 
in efficiency of 7r° invariant mass window requirement be¬ 
tween data and signal MC simulation. It is determined 
to be 1.4% for the rf —y uj'y mode. Since the 7 r° kinemat¬ 
ics in the rf —y uje + e~ decay is similar to the rf —y 077 
mode, the same value is taken as the uncertainty from 
this source for both decay modes. 

The contributions of systematic uncertainties studied 
above and the uncertainties from the branching fractions 
(J/ip —> 777 ' and uj —y 7 T" t " 7 r _ 7 r°) and the number of J/ip 
events are summarized in Table II, where the total sys¬ 
tematic uncertainty is obtained by adding the individual 
contributions in quadrature, assuming all sources to be 
independent. 


V. RESULTS 

The signal yields and detection efficiencies used to cal¬ 
culate the branching fractions and the corresponding re¬ 
sults are listed in Table. I. Using the PDG world averages 
of B(J/ip —y 777 ') and B(uj — y [ 1 ], the branching 

fractions of rj' —y uj'y and rf —y uje + e~ are determined to 
be B(rf —y uj'y) = (2.55 ± 0.03(stat)±0.16(syst)) x 10 “ 2 
and B(if -A uoe + e~) = (1.97 ± 0.34(stat)±0.17(syst)) x 

10~ 4 , respectively. The ratio ^ is then deter¬ 

mined to be (7.71 ± 1.34(stat)±0.54(syst)) x 10 -3 , where 
several systematic uncertainties cancel, e.g., the uncer¬ 
tainties associated with the charged pions (MDC track¬ 
ing), photon detection efficiency, branching fractions of 
J/ip —y 777 ' and uj —y Tr + TT~n° and the 7 r° mass window 
requirement. 


VI. SUMMARY 

With a sample of 1.31 billion J/ip events collected 
with the BESIII detector, we have analyzed the de¬ 
cays rf —y 077 and rf —y uje + e~ via J/ip —> 777 '. For 









the first time, the decay of rf — > we + e" is observed 
with a statistical significance of 8a, and its branch¬ 
ing fraction is measured to be B(r{ —> ue + e~) = 
(1.97 ± 0.34(stat)±0.17(syst)) x 10~ 1 2 3 4 5 , which is consis¬ 
tent with theoretical prediction, 2.0xl0 -4 [2]. The 
branching fraction of v[ —>- wy is determined to be 
B{rf —> wy) = (2.55 ± 0.03(stat)±0.16(syst)) x 10~ 2 , 
which is in good agreement with the world average value 
in Ref. [1] and the most precise measurement to date. 
In addition, the ratio ^ is determined to be 

(7.71 ± 1.34(stat)±0.54(syst)) x 10" 3 . 
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